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We show that the event excess observed by the IceCube collaboration at TeV-PeV energies, 
usually interpreted as evidence for astrophysical neutrinos, can be explained alternatively by 
the scattering of highly boosted dark matter particles. Specifically, we consider a scenario 
where a ^ 4 PeV scalar dark matter particle (f can decay to a much lighter dark fermion y, 
which in turn scatters off nuclei in the IceCube detector. Besides these events, which are ex¬ 
clusively shower-like, the model also predicts a secondary population of events at O(lOOTeV) 
originating from the 3-body decay (f —> xX^y where o is a pseudoscalar which mediates dark 
matter-Standard Model interactions and whose decay products include neutrinos. This sec¬ 
ondary population also includes track-like events, and both populations together provide 
an excellent ht to the IceCube data. We then argue that a relic abundance of light Dark 
Matter particles %, which may constitute a subdominant component of the Dark Matter in 
the Universe, can have exactly the right properties to explain the observed excess in GeV 
gamma rays from the galactic center region. Our boosted Dark Matter scenario also predicts 
fluxes of 0(10) TeV positrons and O(lOOTeV) photons from 3-body cascade decays of the 
heavy Dark Matter particle (p, and we show how these can be used to constrain parts of 
the viable parameter space of the model. Direct detection limits are weak due to the pseu¬ 
doscalar couplings of X- Accelerator constraints on the pseudoscalar mediator a lead to the 
conclusion that the preferred mass of a is > 10 GeV and that large coupling to b quarks but 
suppressed or vanishing coupling to leptons are preferred. 
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I. INTRODUCTION 

The IceCube experiment at the South Pole has recently made international headlines by dis¬ 
covering an excess of events in the energy range from 30 TeV to 2 PeV m- These events are 
usually interpreted as evidence for a flux of astrophysical neutrinos with a power-law spectrum 
~ E~‘^, originating from the production and subsequent decay of charged pions, kaons, muons and 
neutrons produced in collisions of ultra-high energy charged cosmic rays with protons or photons 
in astrophysical sources. 

Despite the plausibility of this explanation, there are other possibilities. For instance, ultra- 
high energy neutrinos could be produced in non-standard processes such as the decay IlHII] or 
annihilation Hals] of very heavy DM particles, or in the early decay of ultra-massive long lived 
particles [HES]- 

In this paper, we explore another alternative idea, namely that IceCube may be observing dark 
matter (DM) particles with PeV energy directly (as opposed to observing only neutrinos from their 
annihilation or decay). The idea, which has first been put forward in [16], is the following: a heavy 
O(PeV) DM species 4>, which makes up a substantial fraction of the dark matter in the Universe, 
decays to a much lighter species x- The resulting flux of highly boosted x particles scatters on 
nuclei in the IceCube detector and leads to the observed energy deposits E^ep up to few PeV. An 
upper cutoff on U^ep is naturally provided by the mass of 4>, explaining the absence of events above 
a few PeV. 

The idea of direct detection of boosted DM in large volume terrestrial experiments was intro¬ 
duced in ref. m in the context of new light (0(1) GeV) particles produced in the annihilation of 
0(100) GeV DM particles in the galactic halo. The authors focused on electron recoil signatures in 
Super-Kamiokande [T 8 |, Hyper-Kamiokande |T9| and PINGU [20|. In a subsequent paper | 2 T], also 
the possibility of detecting boosted particles from the annihilation of heavy DM captured in the Sun 
has been considered. Such signals can be enhanced if the heavy DM particles are self-interacting, so 
that their capture rate in the Sun is increased | 22 j . Also a model with “dark nucleosynthesis” could 
lead to mildly boosted dark sector particles emerging from the Sun |23| . The detection of boosted 
DM annihilation products at much lower energies < GeV in direct DM detection experiments is 
discussed in |24j . The recoil energy spectrum and the annual modulation signal in this case are 
very distinct from those expected from scattering of ordinary non-relativistic DM. Compared to 
these previous works which focus on boosted DM with energies of < 100 GeV, we study signals at 
even higher energies up to O(PeV), and we consider not only direct DM searches, but also indirect 
signatures which may be very relevant in our model. 

In the context of the conventional neutrino interpretation of the IceGube events, there is an 
ongoing debate about the neutrino flavor ratios required to explain the data. The generic expec¬ 
tation for neutrino production from pion decay is that the flavor composition of the astrophysical 
flux at the source (S) is {ue : : i^t)s = (1 : 2 : 0 ) 5 . After propagation and oscillation, the final 
flux at Earth (E) would have a composition of {ue ■ r/fi '■ t^t)e ~ (1 : 1 : 1)^. The IceCube events 
are categorized as track events and shower events, where the former are mostly from charged 
current (CC) interactions with nucleons in which the produced high energy muon leaves a track 
in the detector. The shower events are attributed either to neutral current (NC) interactions of 
neutrinos or to charged current interactions of and v-r- In several analyses, the flavor ratios of 
the IceCube events have been studied |25H31j . and while in general, the data appears consistent 
with a (1 : 1 : l)^; flavor ratio, a mild lack of has been found. 

A unique feature of the boosted DM scenario is that only shower events are predicted at PeV 
energy, while at the lower energies, the ratio of track events to shower events is similar to what 
is expected in the canonical interpretation of the data in terms of astrophysical neutrinos. The 
reason our model predicts also track events at low energy is that, in addition to the dominant flux 
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of boosted x particles, also a secondary flux of DM-induced neutrinos is expected. It arises when 
the particle that mediates DM-SM interactions—taken to be a pseudoscalar a here—is directly 
produced as final state radiation in the heavy DM decay, cj) —)■ xya, and subsequently decays to 
SM particles. While the primary contribution to the IceCube data from x scattering peaks at 
PeV energies but drops at lower energies due to the properties of the pseudoscalar interaction, the 
secondary neutrino flux peaks at 0(100 TeV) energies. Thus, our scenario is also able to explain 
not only the observed ratio of shower to track events, but also the mild (though not yet significant) 
deficit of events in the intermediate energy range of few x 100 TeV. 

Note that the IceCube collaboration has recently published a new analysis [32], the results of 
which are given separately for events coming from above, i.e. from the southern sky, and from below, 
i.e. from the northern sky. This analysis features a notable, but not yet statistically significant, 
bump in the event spectrum from the southern sky at around 80 TeV. Since the galactic center is 
located in the southern hemisphere, a decaying DM scenario like ours predicts a larger contribution 
from the southern sky than from the northern sky. Thus, this bump could be potentially interpreted 
as being due to the secondary neutrino flux discussed in the previous paragraph, which peaks at 
around 100 TeV. 

In addition to the new window to the high energy Universe opened by IceCube, also observations 
at lower energies ~ GeV have caused a stir recently. Namely, an excess of gamma rays from 
the vicinity of the galactic center was found in Fermi-LAT data, which could be explained by 
DM annihilation |33H35j . A good fit to the Fermi-LAT data is obtained for instance for a 30- 
40 GeV DM particle annihilating to bb with a thermally averaged cross-section of about (crurei) ~ 


10“^° cm'^ sec“^, similar to the annihilation cross section expected for a thermal relic. In our 
scenario, a subdominant primordial population of the light DM species x can naturally provide 
such signal by annihilation through s-channel exchange of the pseudoscalar mediator a. We will 
demonstrate that there is a viable region of parameter space which can explain the Fermi-LAT 
gamma ray signal and the IceGube signal simultaneously. 

In the following, we first introduce our toy model of boosted DM in sec. |T|] and then discuss 
the expected IceCube signals in sec. Ill In particnlar, we show which regions of parameter space 
could explain the recently observed high-energy events. In sec. IV we review mechanisms for 
explaining the observed DM relic density [36] in the boosted DM scenario, and in sec.j^we discuss 
the possibility that the galactic center gamma ray excess is explained by XX annihilation along 
with the IceCube PeV events. We then discuss other constraints on the model in sec. VL in 
particular limits from measurements of the cosmic positron and electron spectrum |37H42| . from 
isotropic diffuse gamma rays [431H3] , from direct detection experiments and from searches for the 


pseudoscalar mediator a in flavor physics experiments and at high energy colliders. We summarize 
and conclude in sec. Em 


II. THE FRAMEWORK 


While most of the qualitative results of this paper apply to any PeV-scale boosted DM model, 
we consider as a specific example a toy model featuring a dark sector that contains two DM 
particles: a heavy real scalar (p with mass ~ O(PeV) and a light Dirac fermion x with mass 
~ 0(10) GeV. We denote the relic abundance of cp by /^Ddm and the relic abundance of x by 
/x^DM, where Ddm — 0.258 is the total dark matter density in the Universe [45]. We will discuss 
in sec. IV how and could be determined in the early Universe. We assume that there are no 


other dark relics besides (p and X; be- we assume f(j, + fx = 1. The dark sector Lagrangian reads 


= ^{d'"(t^){dii(p) - ^ml(p‘^ + ix0X - m^xx - y0x<?^TX • 


( 1 ) 
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TTla 


9x 

’^4'/f<t> 

{aVrel)hb 

fx 

BR3(0—>'XX'3^) Comment 

[GeV] [PeV] [GeV] 



[10^5 s] 

[10-26 cm^/s] 



BP 1 12 4.5 30 

0.86 

0.396 

3.6 

2.8 

0.6 

0.022 Vector-like model only 

BP 2 80 3.9 30 

1.51 

0.462 

1.8 

18 

0.33 

0.026 


Table I. Summary of our two benchmark points (BP), both of which can explain the IceCube event excess and 
the galactic center gamma ray excess. In both models, the pseudoscalar a is assumed to couple dominantly 
to b quarks. We also give the calculated values of the velocity-averaged x annihilation cross section (crurei)^^ 
(relevant for the galactic center gamma ray excess), of the fractional abundance of the light DM species 
= l — and of the branching ratio for the radiative decay (/) —>■ xxo. Note that benchmark point 1 can 
be realized only in the Vector-like quark model since in the MSSM-like and Flipped scenarios, laboratory 
constraints on gy^ are too strong (see sec. VI). 


Here, the coupling constant determines the (j) —)• XX decay rate. We assume to be tiny, 
so that the lifetime of cj) is significantly longer than the age of the Universe. One possible way of 
explaining the smallness of y^p^ could be to envision cj) as a composite particle made up of superheavy 
constituents and held together by a new confining gauge interaction. When this new gauge 
symmetry is broken by a tiny amount, a correspondingly small mixing between the and x could 
be generated. Note that we do not include quartic couplings of (j) or Higgs portal couplings in 
eq. ([^ since these interactions will not be relevant to our phenomenological discussion. A possibly 
problematic term could be an operator of the form but we assume the mechanism that 

suppresses yif,^ also forbids or suppresses this operator. 

The light DM x interacts with SM particles through a pseudoscalar mediator a |461449j . This 
pseudoscalar couples to light DM and Standard Model fermions through the Lagrangian 


Tint = igxaXXbX + i'^9Yf 
f 




( 2 ) 


where and gy^ are real couplings of a to light DM x and to Standard Model fermions /, 
respectively, m/ are the SM fermion masses and v ~ 246 GeV is the vacuum expectation value 
(vev) of the SM Higgs field. While generically all gy^ are free parameters, we will specifically 
consider natural scenarios in which the gy^ are generation-independent. 

Throughout most of the paper, we will consider two benchmark points in the parameter space 
of the model, defined in table The heavy DM mass and lifetime r,^, the light DM mass niy, 
and the couplings g^ and gy^ are chosen such that both the IceCube excess of high energy events 
as well as the galactic center gamma ray excess are explained. We assume the mass of a to satisfy 
nia > 10 GeV since constraints are weak in this case (see sec. VI), thus allowing large couplings 
gy^ to fermions. This is important for the model to fit the IceCube data and is also interesting 
because it allows for a detectable indirect signal from the annihilation of non-relativistic relic x 
particles. 

Since the coupling of the pseudoscalar a to SM fermions in eq. (§ should be considered as 
an effective operator after the spontaneous breaking of electroweak symmetry, we need to discuss 
possible ultraviolet completions for such an operator. We consider here three interesting models 
which can provide such a coupling. 

MSSM-like model. In the first model, the pseudoscalar a mixes with an extended Higgs 
sector, for example with the pseudoscalar A® in a type-H Two Higgs Doublet Model (2HDM), by 
a term of the form iaH\H 2 -|- h.c. [50]. In this case, the Higgs couplings to quarks and leptons are 
the same as in the Minimal Supersymmetric Standard Model (MSSM). We therefore denote this 
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model as MSSM-like. The relations for the couplings between effective operator model and the 
complete renormalizable model read m 

9Yd = 9Ye = - tan /3 sin 6/V2 (3) 

9Yu = - cot /3 sin e/V2, (4) 

where tan/3 = U2/U1 is the ratio of the two Higgs vevs and sin0 is the mixing angle between the 
pseudoscalar a and the boson of the 2HDM. Qy^ and qy^ are the generation-independent 
normalization factors of the Yukawa-like couplings for down-type quarks, leptons and up-type 
quarks, respectively. Since the pseudoscalar a couples to SM fermions only through its mixing 
with A^, all of these couplings are suppressed by sin0. As mentioned in the Introduction, we are 
interested in particular in scenarios with large coupling between the pseudoscalar a and bottom 
quarks to optimally fit the galactic center gamma ray excess. This requires large tan /3 to lift up 
the coupling to down-type quarks. Already at this stage, we can see that the MSSM-like model 
will be constrained by experiments sensitive to anomalous couplings of the charged leptons (which 
are also tan/3-enhanced) and by searches for an extended Higgs sector. As we will see in sec. |VI[ 
these constraints lead to the conclusion that the IceCube events and the Fermi gamma ray excess 
can be simultaneously explained in the MSSM-like model only when the pseudoscalar a is heavy 
{rua > mh/2). 

Flipped model. The second model, which we call Flipped is a flipped Two Higgs Doublet 
Model [5TH56] . This means that one Higgs doublet couples to up quarks and leptons, while the 
other couples to down quarks. The difference between this model and the MSSM-like model is that 
the coupling to leptons in the Flipped scenario is proportional to cot j3 rather than tan /3 and is 
thus suppressed rather than enhanced in the large tan/3 region. Therefore, limits from the lepton 
sector will be significantly weaker. The couplings to up-type quarks and down-type quarks are the 
same as in the MSSM-like model. 

Vector-like quark model. The third model has no extended Higgs sector, and the pseu¬ 
doscalar mediator a does not directly couple to SM quarks. Instead, it couples to new, heavy 
vector-like quarks, which in turn mix with the SM quarks m- Since a has no couplings to leptons 
in this model and since there is no extended Higgs sector, we expect constraints to be weaker than 
in the other two scenarios. However, the mass of the heavy vector-like quark should be large to 
avoid LHC limits. 


III. BOOSTED DARK MATTER IN ICECUBE 


III.l. Primary signal: scattering of boosted DM particles on nuclei 


Highly boosted y particles from the DM decay process (j) —)> XX can scatter on atomic nuclei 
in the IceCube detector through their coupling to the pseudoscalar mediator a (see fig. (a)). At 
the high energies we are interested in, the scattering is deep inelastic. Phenomenologically, this 
process is very similar to neutral current scattering of neutrinos, hence its characteristic signature 
is a shower-like event topology. The deposited (or visible) energy E^ep in this case is the energy of 
the recoil nucleus or its fragments. 

The total number of shower events from y scattering in a given E^ep bin 13^“] is given 

by [3n] 


^sh,NC ^ j, 


rm^/2 


dEy X 
^ dE^ 


»^max 

■*^dep 


M^^iEdep) 


dE, 


dep 


dep 


ISmjsf 


-^^dcTp^Ey, F/dep) ^ ^dCniEy, F/dep; 


dE, 


dep 


dE, 


dep 


( 5 ) 
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Figure 1. The Feynman diagrams for (a) the scattering of light DM particle x on nucleons and (b) the 3- 
body decay </> —)■ yya, which produces a flux of high energy pseudoscalars whose decay products contribute 
to astrophysical neutrino, gamma ray and positron fluxes. 


Here, T is the observation time, mjv is the nucleon mass, dap(^n)/dE^gp is the differential scattering 
cross section on protons (neutrons). M^^(£'dep) is the effective detector mass of IceCube for 
neutral current scattering as a function of iiidep- Details on how we estimate M^^(ii'dep) from 
the effective detector mass as a function of incoming neutrino energy, published by the 

IceCube collaboration [1] are given in appendix |A[ Our estimate of M^^(£'dep) is in agreement with 
the results from ref. |2S|, which found that M'^^E^gp) is universal for NC and CC interactions. 

The flux of boosted light DM particles y has a galactic component and an extragalactic 
component 


dEy 


d^ 


GC 


dE^ 


+ 


dE^ 


The galactic contribution is given by 


dcbGC r 

= / dn. 
dE^ I 


1 


diVv 


/ ds/9halo(r(s,'!/’)) 


= 2.1 X 10“^° cm“2 sec“^ x 


1026 


sec 




1 PeV 


\dEx )’ 


( 6 ) 

( 7 ) 


Here, and are the mass and lifetime of the heavy DM particle (j), respectively, phaio is the DM 
density distribution in the Milky Way, r(s, ip) is the position vector relative to the origin at the 
galactic center, s is the distance along the line of sight and ip is its angular direction. We integrate 
the flux over the solid angle and integrate along the line of the sight s. The energy spectrum of 
boosted X particles is simply dN^/dEy- = 6{E^ — 771^/2). The spectrum of antiparticles, dN^^/dE^ 
is the same. The extragalactic contribution to the flux of x particles is |58j 


dEx 


PpMPc 



1 dN^ 
H{z) JEx 


[{l + z)E^]. 


( 8 ) 


In this expression, H{z) ~ v^DA~+~D))((]r+^j3 is the Hubble expansion rate as a function of 

redshift 2;. It depends on the Hubble constant Hq = i^(0), the dark energy density Da ~ 0.692 and 
the matter density Dm ~ 0.308. The cold dark matter density Ddm is 0.258, and the critical density 
of the Universe pc is given by pc — 4.9 x 10“® GeV/cm^ |l5]. Note that we do not account here 
for attenuation of the x flux due to scattering on the interstellar and intergalactic medium. This 
attenuation is already small for neutrinos [58], and the x scattering cross section is even smaller 
than the neutrino scattering cross section. 
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The differential cross section for x scattering on a proton p (neutron n) of mass mj\/ is 

= V ^_ 4sxmN _ ^ fP(^)(x) X - V lAf P 19) 

dx dE(igp ^ 327rs (s — m% — x’^m^rY — Ax‘^m\m% 4 ^ 

q ^ ^ ' NX spins 


with 


i E 

spins 


2glg\ml{Qy 

v‘^{Q‘^ + mYY ’ 


( 10 ) 


where x is the Bjorken scale variable, s = + x'^rnj^ + 2xmi^E^ is the center of mass energy, and 

QY = 2xmjvT'dep is the momentum transfer in the scattering. E^ is the energy of the incoming 
particle y and the nucleon is assumed to be at rest initially, illdep is the energy transferred to the 
hadronic system in the lab frame during the scattering. We are interested in events with a large 
deposited energy E^^p > 10 TeV in this analysis due to the IceCube energy threshold. 

The factor fq^^\x) in eq. (|^ is the parton distribution function (PDF) for protons (neutrons) 
and quark flavor q. We use the PDFs from NNPDF3.0 [59] , which are valid in the range x G Jl0“®, 1] 
and G [2 GeV^, 10® GeV^] and contain the most recent deep inelastic scattering datajj In the 
calculation, we set the PDFs equal to 0 when is smaller than 2 GeV^. Because the cross section 
is proportional to 1/(Q^ + 'ma)‘^i it becomes large when is small. Our cutoff at low would 
therefore affect the results for rua < GeV. In the following, however, we focus on the mass range 
nia > 10 GeV, and we have checked that in this case the contribution of the <2 GeV^ region 
to the cross section is negligible. If one is interested in an extremely light t-channel mediator with 
<C 2GeV^, then the <C 2GeV^ and re <C 1 region, corresponding to exchange of a nearly 
on-shell a, is important. In this region, the PDF description breaks down and one should instead 
calculate the cross section for a* absorption by protons along the lines of the equivalent photon 
approximation in deep-inelastic scattering of electrons on protons. We have used the central values 
of the PDFs, but have checked that varying them within the error band changes the total cross 
section by only 0(10%). The impact on the energy dependence of the differential cross section is 
also negligible. 


III.2. Secondary signal: neutrino flux from 3-body decays of heavy DM 

As mentioned in the Introduction, the boosted DM scenario predicts not only a population of 
high energy events from the scattering of boosted x particles, but also a contribution at lower 
energy from neutrinos produced in the 3-body decay cj) —)> yya (see fig.[^(b)), followed for instance 
by a —7- bb. Since boosted DM can only explain the IceGube events if y particles can scatter on 
nucleons, a mediator particle like a is always needed and the existence of the 3-body decay process 
is thus very generic. Making a heavy does not significantly influence the 3-body decay rate unless 
rUa becomes comparable to m^. The differential decay width of the 3-body decay is, in the limit 
rua —)• 0 and at leading order in niy’, 

^ XXa) _ alvl^Eg - 2m^Ea 

dEa Ibvr^m^ \ 



^ At a; very close to 1, the NNPDF3.0 PDFs are not smooth. Even thongh the large x region is not important for 
onr results, we do not nse NNPDF at a: > 0.1, but use CTEQ5 [SD] PDFs instead. 
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Figure 2. The energy distribution of pseudoscalar particles a produced in the 3-body decay (p —> XX®- 
The parameter values have been fixed at our benchmark values = 4.5 (3.9) PeV, = 30 GeV, 
Too = 12 GeV(80 GeV) for the solid blue (dashed red) lines. 


where Ea is the energy of a in the rest frame of cp. The branching ratio is 

r3(</> -t xxa) 


BR3((?i -t xxa) = 


^3i<P xxa) + T2i(p^ XX) 

r3(</> -t xxa) 
r2(</>^ xx) 


where 


Stt V 


( 12 ) 

(13) 

(14) 


is the rate of the dominant 2-body decay. In the second line of eq. (13), we have assumed that 
is small so that the 3-body decay width is much smaller than the 2-body decay width. We can 
see from Table |I] that this assumption is satisfied at our benchmark points. We plot the energy 
spectrum of a particles from 3-body decay of cp in fig. 

The decay of a to light quarks or b quarks produces neutrinos after parton showering, hadroniza- 
tion and hadron decay. We take the spectra of the secondary neutrinos from each a decay in the a 
rest frame from [61j and boost them into the laboratory frame by folding with the Ea distribution 
from fig. [62]. Multiplying by BR3((/) —)> xxa) gives us the number dNp/dEy of neutrinos per 
energy interval dE^ per cp decay. The flux of secondary neutrinos is then obtained from equations 
very similar to eqs. 0 and ([^ by simply replacing the factor dN^^/dE^ by dN^/dEy. The strength 
of the indirect signal is proportional to g‘^f(i)/T(j, once the masses and rua are fixed. In 

principle, one might also include a factor of the form exp[—Abs(iii^, z)] in the expression for the 
extragalactic flux to account for the absorption of neutrinos in interactions with the cosmological 
relic neutrino background and with the intergalactic medium |58j . However, these effects are neg¬ 
ligible in our analysis and we therefore do not include such an attenuation factor. Moreover, the 
high energy neutrino flux reaching the detector from below is affected by neutrino interactions dur¬ 
ing passage through the Earth. In particular, at energies above ~ 100 TeV, the neutrino-nucleon 
interaction cross section is so large that the Earth can attenuate the neutrino flux. On the other 
hand, electron and muon neutrinos can be regenerated in the decay of tau leptons produced in I'r 
CC interactions. The net effect of both absorption and regeneration is a reduction of the neutrino 
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Figure 3. The galactic and extragalactic neutrino fluxes from 3-body decay of heavy DM, (j) —> xyo, 
followed by a —>■ bb. We have added up the neutrino and antineutrino fluxes and have also summed over 
neutrino flavors. The horizontal dashed line shows the generic flux expected from astrophysical sources, E~'^, 
normalized such that optimum agreement with the IceCube data is achieved [3] . The model parameters are 
set to the benchmark values given in the plot. 


flux by about 15% at neutrino energies ~ 100 TeV [29], and we therefore neglect this small effect 
in our calculation. 

We plot the expected contributions to the nentrino flux from galactic and extragalactic (j) —)• 
XX + (a —t bb) decays in fig. [^for onr two benchmark points. Since the nentrinos originate mostly 
from meson decays after hadronization of the b quarks, their flavor ratio after propagation is 
naturally (1 : 1 : 1)^;. Therefore, we have summed the different flavors, as well as the neutrino and 
antineutrino fluxes, in fig. We see that the secondary neutrinos are softer by about one order 
of magnitnde compared to the boosted DM particles x- The extragalactic flnx is in general softer 
than the galactic one due to redshift. 

Note that, besides the secondary nentrino flnx, there is also a population of boosted DM events 
from scattering of the x particles prodnced in 3-body decay 4> —)• xx®- We neglect these events 
for the following reasons: first, the 3-body branching ratio is almost two orders of magnitude 
smaller than the 2-body branching ratio. Second, the spectrum of x particles from 3-body decays 
is softer than the one from 2-body decays and would therefore contribute only in a regime with 
larger expected backgronnds. Third, a 3-body decay (p —)■ XX® produces only two x particles, 
but typically more than two neutrinos m- Thus, the flux of x particles from 3-body decay is 
subdominant compared to the secondary neutrino flux. Fourth, the x scattering cross section on 
nncleons is usually smaller than the neutrino charged current cross section. 


III.3. Fitting procedure 

To determine the preferred parameter regions for the boosted DM scenario, we use the log 
likelihood ratio (LLR) method. The LLR is defined as follows: 

Sy^xU 9lU \ 

T(j> T(j) J 



LLR 
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= log 


( Max 

a:G[—00,00] L 


/Gauss(x) n* /Poisson ( Si ( 


+ Bi + X ^Bi 


Oi 


\ 


Max 

a:'E[—00,00] 


/Gauss(^0 Oi /Poisson“1“ 


(15) 


Here, /poisson(Ai|n) = /i"-e ^/n! is the Poisson likelihood function and Si{m<f„ gy^g'^f^/Tcf), t^) , 

Bi and Oi are the predicted signal event rate, the predicted background event rate, and the observed 
event rate in the z-th energy bin, respectively. ABi is the Icr error on the background prediction. 
When the nuisance parameter x is 1 (—1), the error x ABi{x) describes the upper (lower) limits 
of the error band, and when x = 0 the background takes its central value. The term fcanssix) 
corresponds to a normal distribution in x and is the Gaussian pull term for the nuisance parameter 
X. By using only one nuisance parameter, we effectively assume that the background uncertainty 
is correlated between bins. 


III.4. Results 


m 


We show the results of our fit in fig. and compare the best fit points to the IceCube data 
fig-i For the mediator mass nia = 12 GeV (80 GeV), the three panels of fig. |^give the best 
fit points (black (red) “+” signs) and preferred parameter regions (black unshaded contours (red 
shaded contours)) at 1, 2, Scr confidence level. For = 80 GeV, the best fit point, marked by a 
red “x” sign, corresponds to one of our benchmark points from tablewhile for rua = 12 GeV, the 
benchmark point (indicated by the black “x” sign) is slightly shifted compared to the best fit in 
order to be consistent also with the galactic center excess and with all constraints. The larger value 
of nia is particularly interesting for the MSSM-like and Flipped models, where it helps to evade 
important constraints from Bg —)• decays and from h ^ aa decays, (see sec. VI.4). Note that 


we 


parameterize the parameter space in hg. in terms of three parameters: the heavy DM mass 
the combination gY^g‘^f 4 >/T(j, of the a coupling constants, the cosmological abundance of the 
heavy DM particle (j) and its lifetime to which the x scattering rate is proportional; and the ratio 
9xf<t>/'’’(!> to which the interaction rate of secondary neutrinos is proportional. In the upper left hand 


plot, we also show constraints from the diffuse 7 ray flux (see sec. VI.2) as thick black (red) lines. 


We always fix the mass of the light DM particle at = 30 GeV, as motivated by the galactic 
center gamma ray excess, see sec. |Vj As expected, the best fit point is always around ~ 4 PeV 
due to the lack of IceGube events above 2 PeV. In hg. we compare the IceGube data from ref. [3] 
to our predictions at the benchmark points. We also show the individual contributions to the 
spectrum separately: the atmospheric (“ATM”) neutrino background (red dotted), the galactic 
(brown dashed) and extragalactic (black dot-dashed) huxes of boosted x particles, and the hux of 
secondary neutrinos from (j) —)• xx + (a —^ bb) decay (purple dashed). 

We see that both the galactic and extragalactic x fluxes contribute at PeV energies, with the 
latter being somewhat softer due to redshift. Actually, the integrated fluxes of the two components 
are comparable, but since the scattering cross-section is higher when the energy of the incoming x 
particle is larger, the softer component is subleading experimentally. Below 1 PeV, the boosted DM 
event rates drop because of the dependence of the scattering matrix element, eq. (10). In their 
place, the secondary neutrino flux takes over below ~ 500 TeV, so that a good fit to the IceGube 
data is obtained at all energies. Note that the normalization of the secondary neutrino flux is set 
by the parameter combination g'^f(j,/T^ and is thus not directly correlated with the boosted DM 
scattering rate, which is proportional to gy^g’^fcji/Tij,. 

Gomparing our two benchmark values of rUa (shaded vs. unshaded contours in fig. left vs. 
right panel in fig. , we observe that the choice of rua has a small influence on the spectral shape of 
the DM contributions, but its main impact is on the overall rate. Therefore, at larger ma, the best 
















11 





Figure 4. Preferred parameter regions for the boosted DM scenario from our fit to IceCube high energy 
data [3]. The three panels show 2-dimensional projections of the 3-dimensional parameter space spanned 
by the heavy DM mass the product to which the scattering rate of boosted y particles is 

proportional, and the combination to which the flux of secondary neutrinos is proportional. (Here, 

and gy^ are coupling constants, is the cosmological abundance of 4>, and is its lifetime.) Solid black 
unshaded (red dashed shaded) contours show the preferred parameter regions at 1, 2, 3cr for vria = 12 GeV 
(ma = 80 GeV) and the black (red) “-I-” signs indicate the best fit points. At rUa = 80 GeV, the best fit 
point is identical to one of our benchmark points (red “x” sign) from table |Tj while for rria = 12 GeV we 
define our benchmark point (black “x” sign) slightly away from the best fit. This way, both benchmark 
points can also explain the galactic center gamma ray excess and evade all constraints. In the upper left 
hand plot we also show as a thick black (thick red) curve the strongest exclusion limits on the nia = 12 GeV 
(toq = 80 GeV) benchmark model, coming from diffuse 7 ray searches (see sec. VI.2). We use = 30 GeV 
for the mass of the light, boosted, DM particle here, motivated by the galactic center gamma ray excess, 
but note that does not affect the IceCube event rate as long as ^ m^. 
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(a) (b) 

Figure 5. Comparison of IceCube high energy data [3] to the prediction at our two benchmark points (see 
Table |T] We plot the signals from galactic (brown dashed) and extragalactic (black dot-dashed) (j) ^ XX 
decays, as well as the contribution from secondary neutrinos produced in 0 —)■ xx-l-(a —)■ bb) (purple dashed) 
separately. The red dotted lines show the atmospheric neutrino background (“ATM”), the blue bars depict 
the background uncertainty and the solid blue lines show the total expected event rate. We have taken the 
mass of the pseudoscalar mediator TOq to be 12 GeV (80 GeV) in the left panel {right panel). We always 
use = 30 GeV for the mass of the light (boosted) DM particle here, motivated by the galactic center 
gamma ray excess, but note that does not affect the IceCube event rate as long as m^. 


fit value of gy^gl^ffp/Tfj, is significantly larger than at smaller ma- When a is heavy, one either needs 
large gvi^gx coupling to keep the scattering cross section of the boosted DM particle x on nucleons 
unchanged, or the flux of x particles must be enhanced by decreasing the heavy DM lifetime r^. 
Note that the two benchmark models shown in figs and [^explain not only the IceCube data, but 
also the galactic center gamma ray excess (see sec. 0 and are consistent with all constraints (see 
sec. 


VI). 


An interesting aspect of our boosted DM scenario is that a dip in the event spectrum is predicted 
between recoil energies of ~ 400 TeV and 1 PeV. This dip is more pronounced at larger nia, see right 
panel of hg. This is in excellent agreement with the current data, which does not feature any 
events in this energy range. Therefore, if this lack of events should become statistically signihcant 
in the future, the boosted DM scenario wonld provide one possible explanation of it. Another 
interesting aspect of our scenario is that, at low energies, where the flux is dominated by neutrinos, 
the expected flavor ratio is (1 : 1 : 1)^; after propagation for most decay modes of a. Thus the ratio 
of shower and track events is predicted to be the same as for the canonical astrophysical neutrino 
interpretation at T^dep ^ few x 100 TeV. On the other hand, at T^dep ~ 1 PeV, the predicted event 
rate is entirely dominated by the DM contribution, which only provides shower events. This is a 
unique feature of this model and can be tested with future data. 


Let us also remark that a recent IceCube analysis |32] which separates events from the northern 
sky and from the southern sky, exhibits a noticeable, but not yet statistically significant, bump 
at energy deposits around 80 TeV in the southern sky. If this bump should become significant 
in the future, it could be interpreted as being due to a relatively large secondary neutrino flux in 
the boosted DM scenario. Since the galactic center, from where most of these secondary neutrinos 
are expected to come, is located in the southern sky, and because neutrinos from the northern 
hemisphere suffer some attenuation in the Earth, our model could explain why a similar bump is 
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not observed in the northern sky. 

Note that, without the neutrinos from the 3-body decay (j) —)• XX + (a —^ bb), the IceCube fit 
of our boosted DM scenario becomes much worse because the prediction would fall short of the 
observed number of events at energies ~ 100 TeV. This could be avoided if a mediator with scalar 
rather than pseudoscalar couplings to fermions, or a vector boson mediator is considered. In this 
case, the boosted DM scattering cross section would not be proportional to (Q^)^, and scattering 
of X particles could explain the IceCube event excess across the spectrum. However, as we will 
argue in sec. VI.3, direct detection constraints in this case may be prohibitive. Ways to avoid these 
constraints include models with inelastic DM scattering or with a very small < 3 GeV, below 
the direct detection threshold. The second possibility would preclude a simultaneous explanation 
of the IceCube events and the galactic center gamma ray excess. 


Let us finally discuss the morphology of the IceCube signal from boosted DM. While the extra- 
galactic flux /{dEy^ dXl^) is isotropic, the galactic component d^^^/{dEy^ dQ^) peaks in the 
galactic center region. (Here V' denotes the direction of sight.) The angular resolution in IceCube 
is about I0°-20° for shower events [1]. With this resolution and more statistics, a morphology 
study of the high energy events would provide an important consistency check of the boosted DM 
hypothesis. 


IV. DARK MATTER RELIC DENSITY 


An important problem of the boosted DM scenario which we have not addressed yet is how a 
particle with a mass of order PeV can account for the observed DM density in the Universe. For 
instance, thermal freeze-out is not a possibility at masses above few hundred TeV due to unitarity 
constraints [63]. A long-lived dark matter particle with a mass of O(PeV) can nevertheless have 
the correct abundance in the Universe |64H68j . 

Non-thermal production mechanisms for PeV DM include jnZj: (I) production in cascade decays 
of the inflaton. In this mechanism, the DM abundance depends on the number density of inflatons 
and on the branching ratio of inflaton decay to DM. (2) production through inelastic scattering 
between high energy particles from inflaton decay and the hot plasma. When high-energy daughter 
particles scatter on the thermalized plasma, DM can be produced until the daughter particles’ en¬ 
ergy become less than Et^ = (dT). (3) For low reheating temperature, DM could be thermally 

produced with the correct relic abundance even when the maximum temperature of the Universe 
during reheating, Tmax, is larger than as long as the reheating temperature (defined as the 
temperature at which the inflaton energy density equals the radiation energy density) is smaller 
than m^. The reason is that the continuing decays of the inflaton produce entropy after DM 
freeze-out, diluting the DM abundance. The authors of ref. m show that these mechanisms can 
account for the abundance of DM with C)(PeV) mass. Mechanism (2) can achieve this even if the 
inflaton does not decay to DM and is thus highly model independent. PeV DM (j) produced through 
this mechanism can for instance account for the observed abundance of DM in the Universe if the 
reheating temperature of order 10 GeV and the mass of inflaton is of order 10^® GeV. |67j. 

In addition to the non-thermally produced relic abundance of heavy DM particles (f>, there 
could also be a thermally produced population of the light DM species y if the thermally averaged 
cross section ((Tt^rei) for xX annihilation through s-channel exchange of the mediator a is not too 
large. This is naturally realized in our scenario, {crv^ei) receives contributions from two classes 
of processes, shown in fig. annihilation to ff and, if rua < also annihilation to aa. The 
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Figure 6. The Feynman diagrams for annihilation of the light DM particle y into (a) SM fermions and (b) 
light pseudoscalar mediator particles a. (The second process is only possible if nia < 'rn^.) 


thermally averaged annihilation cross sections read |48) 


-^27r (4m.; -mi)^ + V ' 


(^^rel)aa 


2tt (4 m 2 — + : 

g^rriy^ {m^ - 6T 

247r {2rri^ — 


2 

x ’ 


(16) 

(17) 


where ruf are the SM fermion masses, the sum runs over all SM fermions /, is the total decay 
width of a, the color factor nI is 3 if / is a quark and 1 if / is a lepton, and T is the temperature. 
The thermally averaged cross section for annihilation to leptons is completely analogous to eq. (16) 
except for the color factor. Note that eqs. (16) and © are approximate results, with only the 
leading terms in the relative velocity Urei kept. The proportionality to T in eq. © arises because 
the process XX clcl is p-wave suppressed. When evaluating for calculating the relic 

density of y, we set T to its typical value at freeze-out: Tp ~ m^/20 [69]. Due to the temperature 
dependence, annihilation to aa can be important in determining the thermal relic abundance of 
X, but does not lead to observable indirect signals today, where the relic population of x is non- 
relativistic. xX f f i th® other hand, is an s-wave process and is therefore relevant both today 
and in the early Universe. 


At our first benchmark point from table |T] (tUq = 12 GeV), it is indeed the interplay of the 
annihilation processes xX QQ and xX bb that sets the relic density of X: fx — O-®- 
second benchmark point {nia = 80 GeV), annihilation to aa is kinematically forbidden at freeze-out, 
therefore xX bb accounts for the relic density ~ 0.33 alone. 


In fact, the thermal production of x has some subtlety to it if the abundance of the heavy 
species (j) is explained by a low reheating temperature Trh- The freeze-out temperature Tp of x 
is of order Tp ~ my^f20 ~ 1.5 GeV at our benchmark points. If Trh ^ Tp, the relic abundance 
of X will be smaller than predicted from the naive estimate for Dirac fermions, ~ 6 x 

10^^ cm^ sec“^/ {av^ei}- If Trh Tp, the thermal production of x is not affected. This is possible 
with a 10^^ GeV inflaton field with Trh 10 GeV that could provide the correct relic abundance 
for (j) |67|. For simplicity, we assume in the following that this second case is realized. We moreover 
assume in the following that cj) and x have comparable relic density, and that together they account 
for all the DM in the Universe (i.e. frf, + = 1). 
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V. THE GALACTIC CENTER GAMMA RAY EXCESS 


The fact that the light DM species x in the boosted DM scenario can have a non-negligible 
relic abundance and a relatively large annihilation cross section to SM fermions in the present day 
Universe indicates that there may be interesting indirect signatures, in addition to the primary 
signal from highly boosted x particles from (p decay. 

In particular, the boosted DM scenario can fit the excess of gamma rays which has been observed 
from the direction of the galactic center at energies of few GeV [33H35] . It has been argued that, if 
the dominant DM annihilation channel is XX bb, as in our boosted DM scenario, a 30-40 GeV 
DM particle with {avre\)i,i in the range 1.4-2.0 x 10“^®cm^ sec“^ provides a good fit to the data. 
Since in our scenario the light DM species x constitutes only a fraction of the total DM relic 
density, its annihilation cross section today has to be correspondingly larger by 1//^. 

At our benchmark points from table the predicted annihilation cross sections are {avre\)j,i ~ 
2.8 X 10“^® (18 X 10“^®) cm^/sec. Here the first number stands for the benchmark point with 
rua = 12 GeV, while the second one (in parenthesis) is for the benchmark point with nia = 80 GeV. 
With fy- = 0.6 (0.33) (see sec. IV), and taking into account that we chose my^ ~ 30 GeV at the 
benchmark points, we thus see that the galactic center gamma ray excess could be explained by 
our boosted DM scenario. Note that for the special case rUa ~ 2m;^, this could be achieved even 
for much smaller couplings gy,, and because the annihilation would be resonantly enhanced. 


VI. CONSTRAINTS 


Constraints on the boosted DM scenario arise on the one hand from indirect DM searches 
sensitive to high-energy particles from the 3-body decay cp —)• XX^, followed by decay of the 
mediator a into SM particles including positrons and gamma rays. We will discuss these possibilities 
in secs. VI. 1 and VI.2, respectively. On the other hand, direct DM searches could hope to directly 
observe the relic population of light DM particles X) see sec. VI.3 Finally, the mediator a could 


be directly produced in accelerator experiments, leading to constraints as well 


see sec. 


VI.4). 


VI. 1. Positron flux from 3-body decay 


XXa 


The flux at any given point x in the galaxy is given by [6T 


d^e±iEe,x) 

dEe, 


1 P(x) 

b{Ee,x) 


XXa) 


^dN^^(E^) 

dES ey 


dES 


I{Ee,E^,x) , 
(18) 


where p{x.) gives the DM density distribution in the galaxy, r3((/) —)• xX®) is the 3-body decay 
rate from eq. E^ is the energy at production, and dN^^{E^)/dE^ is the spectrum at 

production for a decay to //. We obtain dN^±{E^)/dEf in analogy to the secondary neutrino 


spectrum discussed in sec. III.2 by folding the spectrum in the a rest frame (taken from |61j i 
with the energy distribution of a particles from (p —)• xx^ (see eq. 0 and fig. [^. The sum in 
eq. (18) runs over all final states of a decay, and BR(a —)■ //) are the corresponding branching 


IS 


ratios. The factor 6(£'e,x) describes energy loss during propagation |6l]. Finally, I{Ee, E^,x) 
the generalized halo function, which can be understood as a Green’s function of the diffusion-loss 
equation, describing the probability for an with initial energy E^ to be detected with energy 
Ee- We take the halo function from ref. m, assuming a Navarro-Frenk-White (NFW) DM density 
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(a) (b) 

Figure 7. The positron flux from (j) —^ decay, where a decays dominantly to bb. The parameters in the 

left panel (right panel) are fixed at = 4.5 PeV (3.9 PeV) for the heavy DM mass, = 30 GeV for the 
light DM mass, and rua = 12 GeV (80 GeV) for the mediator mass. The AMS-02 positron flux data [39] . 
as well as the Fermi-LAT [?D] and H.E.S.S. [HI |32] data for the combined electron plus positron flux are 
plotted as well. 


profile [70] and the MED propagation model m- The dependence of our results on the DM density 
profile is quite small because the dark matter decay rate only depends linearly on the DM density. 
The uncertainty from the propagation model could change our constraints, but we have checked 
that even for the propagation model MAX from ref. mi, the predicted flux is at most a factor of 
2 larger than for the MED model. 

In fig. we have plotted the positron flux at Earth from the (j) —)• yya decay, where a dominantly 
decays into bb. We fix the mass parameters at our benchmark values nifp = 4.5 PeV (3.9 PeV), 
= 30 GeV and rUa = 12 GeV (80 GeV) in the left panel (right panel). Once the masses are 
fixed, d^f,±/dEe depends on the model parameters through the ratio 

The background model for the e"*" flux is taken from refs. [Z3E3], while the background model 
for the combined e"*" + e~ flux is taken as a fitting function from ref. m- We compare to the 
AMS-02 e"*“ flux data |3S] as well as the Fermi-LAT (30] and H.E.S.S. |41l 132] -|- e~ flux data to 

provide a constraint on this decay. Note that when comparing to Fermi-LAT and H.E.S.S. data, 
which includes both e"*" and e~, the signal flux is twice the signal flux. The error bars in the 
H.E.S.S. data do not contain systematic uncertainties, while those in the Fermi-LAT and AMS-02 
data do. By requiring that the signal flux should be outside the la error bar for any of these data 
points, we find constraints on the coupling g^, the relative abundance of the heavy DM f^, and its 
lifetime T^j,: 


^ < 4.4 X 10-26 sec-i 
< 3.5 X 10-26 sec-^ 

We see from table |I| that our two benchmark points easily satisfy these constraints. 

The cosmic electron background is complicated and model dependent. The background model 
from ref. [H] has a lot of parametric freedom regarding in particular the overall normalization, 
which could alleviate the constraints. Our constraints should therefore be considered as very 


for iria = 12 GeV, 
for nia = 80 GeV. 


(19) 
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conservative. Even for the most conservative assumption of zero background, we would still obtain 
a constraint on by requiring that the predicted signal does not significantly overshoot the 

data. The dominant constraint in this case would come from the last two bins of H.E.S.S. data, 
and the constraint would be weaker by a factor of ~ 5 compared to eq. (19). Also including the 
systematic error of the H.E.S.S. data would make the constraint even weaker. 


VI.2. Gamma ray flux from 3-body decay 


XXa 


The secondary gamma ray flux from the decay 4> —)• yya may contribute to gamma ray searches, 
in particular to gamma ray searches in the galactic center region and in measurements of the diffuse 
isotropic gamma ray flux, i.e. the residual flux obtained after subtracting the contribution from 
known astrophysical sources. We focus here on the diffuse flux because we will see that the 
strongest limits are coming from air shower detectors located in the northern hemisphere and thus 
unable to observe the galactic center m- The only exception is a 7 ray search carried out by 
the IceCube collaboration using the IceTop array m- This search, however, is only sensitive at 
energies above 1 PeV, where the secondary 7 ray flux from decay of ~ 4 PeV DM particles is 
already negligible. Moreover, it is worth emphasizing that searching for signals of decaying DM 
in DM-rich, but also foreground-rich, regions like the galactic center is much less promising than 
searching for annihilating DM in these regions. The reason is that the DM decay rate depends 
linearly on the DM density /?(x), while the annihilation rate scales as /o(x)^. 

The procedure for calculating the diffuse gamma ray flux is similar to the one for the secondary 
neutrino fluxes described in sec. Ill and for the fluxes described in sec. VI.1 In particular, we 
can use eqs. ([^ and ([^ after replacing by the 7 energy and the DM spectrum dN^/dEy_ by 
the gamma ray spectrum at production dN^/dE^. Note that dN^/dE^ must be normalized such 
that its integral over E^ gives the average number of photons produced in each o decay, accounting 
for two body decays without photon emission and for three body decays that lead to the radiation 
of photons. We obtain dN^/dE^ by boosting the 7 ray spectra in the a rest frame (taken from |61) ) 
into the lab frame according to the energy spectrum of a particles given by eq. 0 and fig. and 
multiplying by BR 3 ((/> —)> xyo) from eq. (13). Eor the gamma ray flux, also an absorption factor of 
the form 


exp[—Abs(E^, z)] ( 20 ) 

must be included in eq. (j^ to describe the attenuation of extragalactic gamma rays on their way 
from the source to us. We take this factor from ref. m- 

We then obtain the diffuse gamma ray flux conservatively according to the formula m 


^^diffuse 

dE.y 


d^EG , . d^GC 

-h dvr- 

dE^ dE^ dO. 


minimum 


( 21 ) 


Here, d4>Gc/(c?.£'7 ^^^) I minimum denotes the minimum of the differential galactic flux over solid 
angles, which we take to be the flux from the direction opposite to the galactic center ITS]. We have 
checked that using instead the average of the differential flux over a cone with opening angle 90°, 
centered around the direction opposite to the galactic center, would change d^Qc/{dE^ dD)| min im iim 
by 0 ( 20 %). 

We plot the galactic and extragalactic contributions to the diffuse gamma ray flux in fig. We 
see that the contribution from (p decay in the galaxy dominates over the extragalactic flux due to 
the attenuation factor eq. ( 20 ), which suppresses the extragalactic gamma ray flux. 

Note that we neglect the low energy contribution from inverse Compton scattering (ICS) of 
high-energy from the decay of heavy DM p on CMB photons, starlight, and light rescattered 
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(a) (b) 

Figure 8. The diffuse galactic (solid purple) and extragalactic (dashed brown) gamma ray fluxes from 
(j) —> xya decay, followed by a —> bb. The galactic flux is assumed to have in every direction the magnitude 
it has in the direction opposite to the galactic center m, evaluated assuming a Navarro-Frenk-White DM 
density profile cni. We include only prompt gamma rays, neglecting the low energy contribution from 
inverse Compton scattering because we have checked that the limit is dominated by the prompt signal. 
We compare to the Fermi-LAT measurement of the diffuse gamma ray flux from ref. [43) . using foreground 
model C defined in this reference, and to the limits from air shower detectors The model parameters 

are fixed at the values given by our first (second) benchmark point from tablejljin the left panel [right panel). 


on dust. We estimate [sniEl] that the energy spectrum of ICS photons induced by cf decay peaks 
at 1-100 GeV. Following [9], we have then estimated that the energy density in ICS gamma rays 
predicted at our benchmark points is at least one order of magnitude lower than the energy density 
measured by Fermi-LAT at 1-100 GeV UM- Similarly, also the contribution from bremsstrahlung 
of on dust is negligible. 

To set limits on the parameter space of boosted DM, we compare to the diffuse gamma ray 
spectra from Fermi-LAT |43j and to the flux limits from the air shower detectors KASCADE 
GRAPES-3 [78] and GAMMA [79] , see also [Tlj . From fig. we see that the constraint will come 
mostly from the air shower detectors and the last bin of Fermi-LAT data. By requiring that the 
predicted signal is smaller than the limit from the air shower detectors, we obtain the constraints 


< 0.76 X 10 sec ^ for nia = 12 GeV, 

< 1,44 X 10“^® sec“^ for nia = 80 GeV. 

We see that both of our benchmark points from table [I] satisfy these constraint. 


( 22 ) 


VI.3. Direct detection 

In the boosted DM scenario, conventional DM direct detection experiments can only constrain 
the thermally produced population of light DM particles X; iiot the population of heavy DM 
particles cf. The density of (f particles and thus also the flux of boosted x particles from cf decay 
are too small to be observed in these detectors. Therefore our discussion of direct detection will 
focus on the non-relativistic population of the light DM species x- The cross section for x^aucleus 
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scattering is [18] 


da 

dEr 


i^T 9x (Q^)^ 


( 23 ) 


where is the nuclear recoil energy, v is the DM velocity, = 2mTEr 100 MeV^ is the 
4-momentum transfer squared, rriT is the mass of the target nucleus and mjsf is the nucleon mass. 
The quantities E^,’,^ are the pseudoscalar form factors of the target nucleus (see e.g. [82]), and 
the effective nucleon couplings g^, g'jx depend on the gvf (see also ref. [IH]). For our choice 
^ 10 GeV, we have ;§> Q^, so that is negligible in the denominator. The factor (Q^)^ 
in the numerator arises because, in the non-relativistic limit, xi^X ^ '\f^- Direct detection 
constraints are in general very weak in our boosted DM model due to the (Q^)^ suppression unless 
the mediator mass ma is extremely small. The resulting limit on gx9Yf is therefore much weaker 
than the value needed by the thermal relic density |48j . 

Departing for a moment from our toy model with a pseudoscalar mediator, we note that in 
general, boosted DM models with interaction cross sections strong enough to explain the IceCube 
events would also lead to a large signal in direct detection experiments. From a model building 
point of view, there are several ways of circumventing this, other than using a pseudoscalar coupling 
as in our toy model. (1) Construct a model in which the scattering of the light DM particles x 
on nuclei is inelastic [H3]. If the mass splitting 5m between the ground state of x and the excited 
state X* which is produced in the scattering is sufficiently large, it will lead to vanishing event 
rates in direct searches, but will have no influence on boosted DM collisions as long as 5m is small 
compared to the energy of the boosted DM particles. (2) Assume the relic abundance of the light 
DM species is sufficiently low to avoid direct detection limits. This would of course preclude a 
simultaneous explanation of the IceCube events and the galactic center gamma ray excess. (3) 
Choose the light DM mass smaller than ~ 3 GeV, below the energy threshold for direct detection. 
This would also preclude an explanation of the galactic center gamma ray excess. 


VI.4. Constraints from flavor physics experiments and from collider searches 

In the following, we discuss constraints on our boosted DM scenario from experiments at fla¬ 
vor factories and at high energy colliders and indicate for each constraint to which of the three 
renormalizable models from sec. |I^it applies. 

A large number of constraints arises from Kaon and B meson decays |48j . Searches are sensitive 
to the production of the pseudoscalar a in decays of these mesons if a subsequently decays to 
leptons, photons or invisible particles. Since we are considering the case ma > 10 GeV, those 
constraints are, however, significantly weakened by the fact that a would have to be off-shell. 

Bs is the only search channel sensitive to an off-shell pseudoscalar. If we consider a 

renormalizable model for the pseudoscalar a in the framework of a Two Higgs Doublet Model, as 
in the MSSM-like and Flipped models from sec. [ITJ a couples to the SM by mixing with the heavy 
pseudoscalar A^. The mixing angle is denoted by 6. The branching ratio for Bg —?• in the 

MSSM-like model is given in ref. [5D1IH1]. The contribution from a to the amplitude is proportional 
to tan^/3sin^0. The constraint for ma ~ 10 GeV is tan/3 sin 0 = ^j2gYjj^ < 0.4 (0.51) for 
charged Higgs boson masses oi mH± ~ 800 (400) GeV, while the constraint for tUq ~ 80 GeV is 
about tan/3 sin 0 = ^2gY^gY^ ^ 3.8 (4.8) [50]. For the Flipped model, where lepton couplings are 
proportional to cot /3, the amplitude from a exchange is proportional to tan /3 cot /3 sin^ 9 = sin^ 6. 
Therefore, the constraint is sin0 = x/2gY^gY^ ^ 0.4 (0.51) for charged Higgs boson masses of 
mjj± ~ 800 (400) GeV when ma ~ 10 GeV. For ma ~ 80 GeV, there is no constraint on sin0. 
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Note that the Bg —)• constraint does not apply to the Vector-quark model because a does 

not couple to leptons in this model. 

An additional constraint, which is independent of the couplings of the pseudoscalar a to 
fermions, arises from the exotic decay h —)• aa. In the context of the MSSM-Uke and Flipped 
models, the branching ratio for this decay is constrained by [50l |85l | 86 ] 

BR(/. 0.02 (^^)‘(^)‘< 0.22, (24) 

where tua is the mass of the heavy pseudoscalar. If ruA — 800 GeV, sin 6 has to be smaller than 
0.02. If nia becomes comparable to m/i/2, the above constraint is weakened, and for nia > mh/2 
it is completely absent. It is also absent in the Vector-quark model. 

We should also consider constraints from the LEP experiments, which have searched for e'’'e“ —)• 
hA^ , where is the pseudoscalar Higgs boson appearing in the MSSM [HZj. While these searches 
exclude A^ masses below 90 GeV, they do not apply to models with an extra pseudoscalar a, like 
the scenarios we are considering here | 88 j . 

If a is heavy enough to decay to XX; ref. m shows that searches for b jets and missing energy 

can provide an excellent constraint on the pseudoscalar a. The dominant processes are gg —)• bba 
(—) (—) 

and b g ^ b a, with a decaying to XX subsequently. The current CMS and ATLAS searches 
[ 8 ^ [90] . which are optimized for final states with two b quarks, lead to the constraint ^ 5 

for nia ~ 100-250 GeV and assuming g^ = gY^\/2mi,jv [57]. If gy. is significantly larger than 
gYf^y/2mbjv, the limit will become somewhat weaker since the probability for radiating an on-shell 
a particle changes m- 

In the intermediate mass region 20-80 GeV, ref. [ 88 ] also discusses the processes gg —)• bba and 
(—) (—) 

bg ^ b a, but considering the subsequent decays a —>■ and a —)• r+r”. By looking for 

these leptonic final states, the high luminosity LHC can be sensitive to gY^ 7 with 100 fb ^ of 
integrated luminosity, assuming that gYf is universal for down type quarks and charged leptons (as 
in the MSSM-like model). Since this assumption is not satisfied in the Flipped model, which has 
suppressed couplings of a to leptons, and in the Vector-quark model, in which a does not couple 
to leptons at tree level, the constraint would be significantly weaker or completely absent in these 
models. 

In the light mass region nia ~ 5.5-14 GeV, CMS has searched for a —>■ g~ in the context of 
the Next-to-Minimal Supersymmetric Standard Model [9T]. The upper limit on the cross section 
for the process pp —?• a —)• g'^g~ is around 2-4 pb. This translates into a constraint of gY^ ~ 2 
in the MSSM-like model, where = gY^ | 88 j . The Flipped and Vector-quark models are not 
restricted by this constraint due to the smallness or complete absence of leptonic couplings of a. 

Let us summarize the most stringent constraints for the three models defined in sec. [H] (see 
also the last column in table [lT| below). For the MSSM-like model, the most stringent limit comes 
from Bg —7- g~^g~- It rules out the MSSM-like model as a UV-completion for our nia = 12 GeV 
benchmark point, while for the nia = 80 GeV benchmark point, it is a viable possibility. For the 
Flipped model, the coupling between leptons and the pseudoscalar a is suppressed once we are in 
the large tan/3 region. But the constraint from h ^ aa still implies that the mixing angle sin0 
between a and the heavy pseudoscalar A^ should be very small. If we require that tan/3 < 50, 
this disfavored also the Flipped model as a UV completion for our nia = 12 GeV benchmark point. 
At nia = 80 GeV, the h ^ aa constraint is absent because the decay is kinematically forbidden. 
For the Vector-quark model, only the perturbativity of the Yukawa couplings involving the heavy 
quarks, together with the LHC limits on their mass, imposes a very weak constraint gY^ ^ 20 m- 
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Table II. Summary of constraints on the boosted DM scenario for two different benchmark values for the mass 
TOa of the pseudoscalar that mediates interactions between the light DM species x and the SM. Since IceCube 
sees both the scattering of highly boosted y particles and secondary neutrinos from ^ yy + (a — >■ bb) 
decay, the experiment constrains two independent combinations of the pseudoscalar couplings to DM 
and b quarks {gYb)j the lifetime of the heavy DM particle, t^, and its fractional abundance in the Universe, 
/ 0 . Note that we always assume here that a couplings to SM fermions other than the b quark are negligible. 
Requiring that the galactic center gamma ray excess can be explained constrains the light DM mass and 
an additional combination of coupling constants. Further constraints come from secondary and 7 rays 
from (j) —>■ XX + (a —> bb) and from laboratory searches for the pseudoscalar mediator a. 


VII. SUMMARY AND CONCLUSION 

In summary, we have discussed the possibility that the high energy event excess observed by 
the IceCube collaboration is explained by the scattering of highly boosted DM particles on atomic 
nuclei in the detector. We have constructed a simple toy model in which a DM particle (j) with a 
mass of order PeV can decay into a much lighter DM species x- The x particles, in turn, interact 
with atomic nuclei through a f-channel mediator a, thus explaining the IceCube signal at PeV 
energies. 

The experimental constraints on this toy model are summarized in table for two different 
benchmark values of the pseudoscalar mass ma. At both benchmark points, we have assumed 
that the mediator a has significant coupling to h quarks, while its couplings to light quarks and to 
leptons are suppressed. This is naturally realized in UV-complete models with either an extended 
Higgs sector or with the introduction of vector-like quarks (see sec. H- The highest energy events 
in IceCube set the scale for the heavy DM mass m^j) and the normalization of the scattering cross 
section. At lower energy, IceCube is sensitive to the secondary neutrino flux from the 3-body decay 
^ XX + {a ^ bb), see fig. B This provides a constraint on the branching ratio for this decay. 
Since the same decay also leads to secondary electron/positron and gamma ray fluxes, and 7 
ray data from AMS-02, Fermi-LAT, HESS and several air shower arrays provide a constraint on 
its branching ratio as well. Moreover, the boosted DM scenario is constrained by searches for the 
new pseudoscalar particle a in flavor physics experiments and at high energy colliders. 

We have shown that, besides explaining a population of high energy events in IceCube, the 
boosted DM scenario can simultaneously also account for the gamma ray excess observed in Fermi- 
LAT data from the direction of the galactic center. This is possible because the light DM species 
X can have a non-negligible thermally produced relic abundance, and can annihilate in the Milky 
Way today. Fermi-LAT data then identifies a preferred range for the light DM mass and its 
couplings to ordinary matter. 
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The boosted DM scenario shares some features with interpretations of the IceCube data in terms 
of DM decay directly to SM particles, including neutrinos. First, the morphology of the signal is 
similar in the two cases, with a mild peak expected in the galactic center region. Moreover, it 
is worth mentioning that the most recent IceCube data [32] provides a mild hint at a bump-like 
feature at ~ 80 TeV from the southern sky. Since this is where the galactic center is located, 
such a bump could be explained by the secondary neutrino flux in the boosted DM scenario. The 
second common feature between boosted DM and more conventional decaying DM explanations 
of the IceCube data is the rapid drop of the signal at energies larger than half of the heavy DM 
mass. With more statistics collected, these features can help to distinguish DM interpretations of 
the IceCube data from an interpretation in terms of isotropic astrophysical neutrino emission. 

A unique feature of the boosted DM scenario is the prediction that, at PeV energies, where the 
IceCube signal is explained by scattering of boosted DM particles, only shower-like events should 
be observed. At lower energies ~ 100 TeV, however, where the secondary neutrino flux from the 
3-body decay cj) —)• XX + (a —>■ bb) contributes, both shower and track-like events are predicted, 
with a ratio very similar to the one expected from astrophysical neutrino sources. Between the two 
populations of events, a mild dip in the energy spectrum is predicted. These features distinguish the 
boosted DM scenario from astrophysical explanations of the IceCube data and from interpretations 
in terms of neutrinos from DM decay. 


ACKNOWLEDGMENTS 

It is a pleasure to thank Carlos Argiielles, Sergio Palomares Ruiz, Hubert Spiesberger and 
Wei Xue for useful discussions. We also gratefully acknowledge discussions with Stefano Morisi, 
who independently had the idea of interpreting the high-energy IceCube events as a DM “direct 
detection” signal. JK and JL are supported by the German Research Foundation (DFG) under 
Grant No. KO 4820/1-1. JK would like to thank CERN for hospitality during the final stages of 
this work. 


Appendix A: Details on the calculation of effective detector mass 

Here, we discuss how we obtain the effective detector mass of IceCube, which 

appears in eq. (§. We use fig. 7 from ref. pQ, which shows the effective detector mass as a function 
of neutrino energy rather than the deposited energy Edep- However, we can exploit the fact that in 
charged current (CC) interactions of electron neutrinos, all the neutrino energy (including both the 
energy transferred to the hadronic system and the energy of the produced electron) is deposited 
in the detector. Electron neutrinos produce shower events very similar to neutrino or DM neutral 
current scattering, hence we can assume 

= (A1) 

We have verified the validity of this assumption by checking that we can use obtained 

this way to reproduce the effective detector mass for NC neutrino interactions according to the 
convolution formula 

Jo \^i') ^-^dep 

where {E ^, Edep) /dEdep is the differential cross section for NC neutrino interaction and 

a^^{Ei^) is the corresponding total cross section [92] . In fig. we compare our result for 
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Figure 9. Effective target mass for neutrino interactions in IceCube as a function of the incoming neutrino 
energy E^, for CC z/g (red solid) and NC (gray solid) interactions [1]. The dashed brown line shows our 
prediction for the effective target mass in the NC case from eq. (A2), which is in excellent agreement with 
the results from [T]. 


with the IceCube data (fig. 7 in [T]), and find excellent agreement. Our results are also in agreement 
with the dedicated fitting result from [29] . 
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